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I. INTRODUCTION 



Stability was originally a physical concept, believed first to have 
been considered by Lagrange in connection with the equilibrium states 
of conservative systems. Equilibrium states are those in which the 
velocities and accelerations are zero, describino a state of rest. 

A system in equilibrium may react in one of two manners when influ- 
enced externally; it may exhibit motions which remain in a small 
reqion about the equilibrium position (which does not preclude the 
possibility that the system return to its state of equilibrium), 
or it may exhibit motions which cause the system to increasingly 
deviate from its equilibrium state, never returning to it. In the 
former case, the equilibrium state is said to be stable, while the 
latter case describes an unstable state of equilibrium. 

The above remarks are, of course, not rigorous and are intended 
solely as an intuitive guide to the meaning of stability of equili- 
brium states of physical systems. Physical systems are almost 
always described by a system of differential equations and in 
studying the behavior of a physical system, we ultimately study the 
behavior of the solutions to the differential equations which 
describe the system. Thus the study of stability of a ohysical 
system is reduced to a consideration of the stability of solutions 
of differential equations, and stability theory is today an integral 
part of that segment of analysis known as the theory of ordinary 
differential equations. 
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stability theory is quite advanced today and there are numerous 
concents and approaches which vary with both the form and order 
of the system of differential equations. For example, the Poincare- 
Bendixson theory is an analytical -toooloqi cal anoroach which is 
particularly useful for second order autonomous systems, while for 
non-autonomous systems, the concept of stabilitv in the sense of 
Lyapunov is utilized and the approach is purely analytical. This 
paper is concerned with a general study of stability in the sense of 
Lyaounov, with particular emphasis on systems of linear differential 
equations though some results for nonlinear systems are considered. 
Section II deals with a review of familiar facts from the theory of 
linear ordinary differential equations. Stability and boundedness for 
linear systems are discussed in section III. Non-homoqeneous systems 
are considered in section IV, and the concents of uniform and restric- 
tive stability are treated in section V. Some results for non-linear 
systems are given in section VI, and concluding remarks follow in 
section VII. . 
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II. PRELIMINARY REMARKS 



We consider the behavior of solutions of the nth order linear 
non-homoqeneous system 
(2.1) x' = A(t)x 

where x = x(t) is an n by 1 column vector and A(t) is an n by n matrix 
with elements a..(t), i , j = 1,...Hj which are real valued functions 

' J 

continuous on the half-line [t„,«>). 

The norm of x, A is defined by 
(2.2i) Mx(t)l| = E|x.(t)l 

(2.2ii) ||A(t)|| = ElA.j(t)| 

and note that both llx|l and ||Al| are functions of t. The usual 
orooerties of norms are satisfied and are stated here without oroof. 
For 1 |x 1 1 they are , 



(2.3i) 

(2.3i i ) 

(2.3i i i ) 
(2.3iv) 

and for | 1 A| | , 



I |x| I ^0 and [ |x| | = 0 iff x(t) = 0, 

llxll - l|y|l| 1 lixtyll illxll + ||y||. 

I |axl 1 = |a| I |x| 1 for a a real number, 

1 l/x(t)dt| I ^ ;| |x(t) I |dt. 



(2.4i ) 11^11^0 and IjAll = 0 iffA(t) = 0, 

(2.4ii) I ||A|| - ||B||| 1 ||A±B|| < ||A|| + ||B|1, 

(2.4iii) IIABII < IIAII ||B||, 

(2.41V) IIAxll 1 llAlMlxll, 

(2.4v) l|ctA|| = |ci| 1 1A| I for cx a real number. 

Canital letters will denote n by n matrices, lower case letters 
n by 1 column vectors, and Greek letters real numbers unless otherwise 
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noted. Primes (') denote differentiation with respect to t; denotes 
the transDOse of A; tr(A) denotes the trace of A; det A denotes the 
determinant of A; and adj A decotes the adjoint of A. Lower subscripts 
denote comoonents of a vector or matrix, for examnle, x^. is the ith 
component of the vector x. 

The questions of existence and uniqueness of solutions to the 
system (2.1) may be dealt with quite briefly. Solutions to (2.1) 
satisfying the initial condition x(t^) = x^ exist and are unique on 
any interval in which A(t) is continuous. Hereafter, continuity of 
A(t) on [Tq,") will be assumed so that existence and uniqueness of 
solutions for systems of this form are always assured. 

Linear homogeneous systems of differential equations are charac- 
terized by the fact that any linear combination of solutions is also 
a solution. The space of all possible solutions to the system (2.1) 
constitutes an n-dimensional vector space. Thus there exists n 
linearly independent solutions, (j)\ (|)^,..., (|)^, which form a basis for 
this vector space of solutions. Any particular solution x(t) of (2.1) 
can then be expressed as a linear combination of elements in this 
basis. Since any n linearly independent solutions constitute a basis, 
it is convenient to pick a particular set and define the fundamental 
matrix X(t) as follows. 

(2.5) X(t) = ^l(t) (i,j = l,...,n), 

i.e., the columns of X(t) are n linearly independent solutions of 
(2.1) and these n solutions (i = l,...,n) are chosen to satisfy the 
initial condition <f>^-(t ) = 6.. (i,j = l,...,n) so that X(t ) = I, 

the identity matrix. 



8 



From (2.5), we note that X(t) satisfies the matrix differential 
equation 

(2.6) X' = A(t)X, 

and that any solution x(t;t^,x^) of (2.1) satisfying the initial 

condition x(t ) = x may be written 
0 0 

(2.7) x(t) = X(t)x^. 

From the formula for the derivative of the determinant of a 
matrix, 

d det X(t) = tr [A(t)] det X(t), 
dt 



we obtain the Jacobi-Liouville formula 

t 

(2.8) det X(t) = det X(t ) exo f tr A(s) ds 

^0 

t 

= exD / tr A(s) ds, 
t 

0 

and note from (2.8) that X(t) is non-sinqular for t ^ t^. 

Any non-singular matrix W(t) which is a solution of (2.6) may 
serve as a basis matrix though if W(t^) /I, we no longer have (2.7) 
but instead 

x(t;to,Xo) = W(t) Wq, 

where w ^ x . 

0 0 

From the formula for the derivative of the inverse of a matrix, 

d [X’’'(t)] = -X’''(t) X'(t) X’''(t), 

-IT 

we obtain the result that if Y(t) = X (t) , then Y(t) satisfies 
the matrix differential equation 

(2.9) Y' = -A(t)^ Y, 
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and since Y(t) is non-singular and Y(t^) = I, we note that Y(t) 
is the fundamental matrix to the adjoint system of (2.1) given by 

(2.10) x' = -A(t)^ X. 

The fundamental matrix X(t) provides an explicit form for the 
solution of a non-homogeneous system of the form 

(2.11) x' = A(t) X + b(t) 
as the following proposition illustrates. 

Proposition 2.1 . The solution x(t;t^,x^) of the system (2.11), 

where A(t) and b(t) are continuous on [t ,«) may be written 

? 

x(t) = X(t) X + / X(t) X“^(s) b(s)ds, 

where X(t) is the fundamental matrix of the homogeneous system 
x' = A(t) X. 

Proof: Direct substitution. 

A similar formula may be obtained for non-linear mixed systems 
of the form 



(2.12) x' = A(t) X + f(t,x). 



Proposition 2.2. Let f(t,x) be continuous for t > t and 
|x|| < c, and A(t) continuous for t ^ t^. Then the solution 



x(t; t ,x ) of (2.12) may be written 

0 0 ^ 



x(t) = X(t) X + / X(t) X" (s) f(s,x(s)) ds, 
° t 



0 

where X(t) is the fundamental matrix of x' = A(t) x. 

Proof: Direct substitution. 

Both Propositions 2.1 and 2.2 simplify in the case A(t) = A, 
a constant matrix and t^ = 0. For then, the matrix X(t) X”^(s) 
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satisfies (2.6) with initial condition 

X(t) X'^s) = I 

at t = s, and the same follows for X(t-s). Hence by uniqueness, 
(2.13) X(t) X‘^(s) = X(t-s), 

for A a constant matrix and X(0) = I. 

Other well known results from the theory of ordinary differential 
equations will be stated when needed. 
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III. LINEAR HOMOGENEOUS SYSTEMS 



Consider a general first order system of ordinary differential 
equations of the form 
(3.1) X' = f(t,x) 

where x(t) is an n by 1 column vector and f(t,x) is an n by 1 
column vector defined and continuous in the region P = {(t,x): 

0^t<«>, I I X 1 1 < . 

We wish to define a concept of stability for this system. One 
definition which immediately comes to mind might be to define the 
system to be stable if it exhibits only finite motions, i.e., all 
solutions remain bounded as t-«». However, such a concent is not 
particularly useful as it distinguishes variational effects only as 
to their being finite of not. For example, an earthquake would be 
stable in this sense. 

What we wish to embody in a stability concent is a means by which 
we can measure the effects of small variations on the system (or a 
particular solution of the system). One criterion we would desire of 
a stable system is that small variations not affect the system in a 
drastic manner. This is the nrimary consideration embodied in the 
concent of Lyanunov stability, which imnoses stringent restrictions 
on the solutions of (3.1), by requiring that solutions which are 
initially close together remain close together for all future time. 
Formally, we define stability as follows. 
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Definition 3.1 . Let x(t;t^,x^) be a solution of (3.1) satisfying 



x(t) exists for t > t^, 

— 0 



(3.2i) 
and 

(3.2ii) 

Then x(t) is said to be stable (in the sense of Lyanunov) if 



(t,x(t)) e P for t ^ t^. 



(3.2iii ) 



there exists a>0 such that every solution 
y(t;t^,yQ) exists on and Ct,y(t)) 



whenever | |x^ - y^ | 



< Ot. 



and 

(3.2iv) 



given e>0 there exists 6 = 6(c,f,x^), 

0<6<a, such that | |x(t) - y(t)| | < £ for t ^ t 



for all solutions y(t) with | |x^ - y^| 



< 6 . 



Strictly speaking, the above defines stability at the right, since 
the behavior as t->«> is considered; replacing the interval by 

the interval (-“,t^] would define stability at the left. Note 
from (3.2iv) that stability is eguivalent to uniform continuity of 
x(t;t ,x ) on [t ,«>) with respect to the initial vector x„. A solution 
which is not stable is said to be unstable. Hereafter, stability will 
always refer to Lyapunov stability at the right as defined above. 

A special case of stability occurs when solutions which are 
initially close approach one another as t^°°. This leads to the 
following definition. 

Definition 3.2. A solution x(t;t ,x ) of (3.1) is said to be 

0 0 

asymototical ly stable on [t^,”) if it is stable and in addition there 
exists 6' > 0 such that | |x(t) - y(t)||-<-0 as t-^ whenever y(t;tQ,y^) 
is a solution with 11^^ - y^M < 6'. 
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Note that by definition, a solution which is asymptotical ly stable 
is stable. If all solutions of the system (3.1) are (asymntotical ly) 
stable, then we say that the system is (asymptotically) stable, while 
if there exists an unstable solution to (3,1), we will call the system 
unstable. 

The following examples illustrate some of these concepts. 

Example 3.1 . Consider the one-dimensional system 

X ' = -X 

on [0,oo). The general solution x(t;0,x^) is given by x(t) = x^ e”^, 
so that 

I |x(t) - y(t)| I = |Xq - .VqI ^ |Xq - y^l 
Conditions (3.2 i,ii, and iii) are satisfied as is (3.2iv) with 
6 = e. Thus the system is stable. 

In addition, 

I |x(t) - y(t) II = |x^ - y^l e'^ ^ 0 as t->~ 

so that the system is asymptotically stable. 

Example 3.2 . Every non-trivial solution of the one-dimensional 
system 



x' = 0 on [0,«>) 

is of the form x(t;0,x ) = x and thus 

0 0 

||>^(t) - y(t)|| - |x^-yj. 

Conditions (3.2i,ii, and iii) are satisfied as is (3.2iv), aaain with 



6 = e . 

However, 

lim I |x(t) - y(t) II = |Xq - y^| > 0, 
t-X» 

so that no solution is asymptotically stable. Thus the system is 
stable. 
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Example 3.3 . Non-trivial solutions of the one-dimensional system 



x' = X on [0,“>) 

are of the form x(t;0,x^) = x^e^ and thus 

I |x(t) - y(t)| I = |Xq - y^|e^ 

which can be made arbitrarily larqe independent of - y^| by taking 
t large enough. Thus the system is unstable. 

Example 3.4 . Cesari [2] shows that the trivial solution z(t) 5 0 
of the one dimensional system 

x' = x^ on [0,®) 

is not stable since every (real) solution x{t;0,x^) is of the form 

x(t) = ^0 

rrag 

and solutions with x > 0 cease to exist at t = (x„)" in violation 

0 0 

of (3.2iii). 

For linear systems of the form 
(3.3) x' = A(t) X + b(t) 

with A(t) and b(t) continuous on [t^,®), (3.2i,ii, and iii) are always 
satisfied so that (3.2iv) becomes the determining factor as to whether 
or not solutions to these systems are stable. For these systems, 
we have the following useful proposition. 

Proposition 3.1 . If the linear system (3.3) has a stable solution 
then all its solutions are stable, i.e., the system is stable. 

Proof: Let x(t;t^,x„) be a stable solution and y(t;t„,y^) 

be any other solution. Let e > 0. 

If z(t;t^,z^) is a solution of (3.3), then x(t) - y(t) + z(t) 
is also a solution since the system is linear. Then since x(t) 
is stable, there exists 6 > 0 such that if 
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l|x(to) - [x(to> - y(to) + z(to)]|| = llyo-"oli 

< 6 

then 

||x(t) - [x(t) - y(t) + z(t)]|| = ||y(t) - z{t)|| 

< e 

which says that y(t) is stable. Thus all solutions are stable. 

Corollary 1 , Pronosition 3.1 remains true if stable is replaced 

by asymptotically stable. 

Proof: The proof is identical to that above. 

We next define the concept of boundedness of solutions to 3.1. 

Definition 3.3 . A solution x(t) of (3.1) is said to be bounded 

if there exists M > 0 such that | | x(t) | | ^ M, t ^ t^, 

i.e., if ||x{t)|| is a bounded function of t on [t^,«) . If all 

solutions to (3.1) are bounded, we say that the system is bounded. 

Stability and boundedness are independent concepts as the 

following examoles illustrate. 

Example 3.5 . All solutions of the one dimensional system 

x' = 2t on [0,“) 

2 

are of the form x(t;t ,x ) = x + t and are stable yet unbounded 

0 0 0 

on [0,~). 

Example 3.6 . Cesari [2] shows that every non-trivial solution of 
the two dimensional system 

X* = y 

y' = - 1/2 X (x^ + (x^ + 4y^)^/^^ 

or equivalently, 

z* ■ = -1/2 z (z^ + (z^ + 4(z')^)^/^ 
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is of the form 



z(t) = c sin (ct+d), c ^ 0 

and is thus bounded; none of these solutions is stable. 

For linear homogeneous systems of the form 
(3.4) x' = A(t) X 

stability and boundedness are equivalent as the following theorem 
illustrates. 

Theorem 3.1 . The system (3.4) is stable iff it is bounded. 

Proof: Assume the system is stable. To show that it is bounded, 

it suffices to show that there exists M > 0 such that ||X(t)|| M 
where X(t) is the fundamental matrix of the system. Since the system 
is stable, in particular, the trivial solution z(t) e 0 is stable. 

Thus if e > 0 there exists 6 > 0 such that every solution x(t;t^,x^) 
with MXqII < 6 satisfies ||x(t)|| < e for t ^ t^. In particular, let 
x(t) be that solution with initial vector x^ given by 




where e^ is the vector whose jth component is given by el = 6... 

J "I J 

Then | | | < 6 so that 

I |x(t)| 1 = I |X(t) Xgl 1 < e, t ^ tg. 

(S i i 

But for this choice of x^, X(t) x^ = ■j- > where 4> is the ith column 

of X(t). Thus 

^ I |(()''(t)| I < e t > t . 

2 ~ ° 

This procedure is valid for i = l,...,n so that 

l|x(t)|| = I ||♦'(t)|| < ^ • M 

i=l ^ 

and hence X(t) is bounded. Thus the system is bounded 
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Conversely, assume the system is bounded so that X(t) ^ M, 
for t ^ t^. Then if x(t;t^,x^) and y(t;t^,y^) are two solutions to 

(3.4) 



I |x(t) - y(t)| I = I |X(t) (Xq - y^)! I 1 M I |x^ - y^ 
so that if I |Xq - y^l I < |- , then | |x(t) - y(t) | | < e for t ^ t^ 
which shows that the system is stable. 

Corollary 1 . The system x' = A(t) x is asymototical ly stable 
iff X(t)->0 as t-»-«>, where X(t) is the fundamental matrix of the system. 

Proof: The reverse implication is obvious. To show the forward 

implication, note that the hypotheses imply that the trivial solution 
is asymptotically stable and the result then follows as in the oroof 
of the first half of Theorem 3.1. 

Before proceeding to consider various results concerning stability 
in linear homogeneous systems, it is convenient to introduce the 
following lemma, known as Gronwall's Lemma, which is an extremely 
useful tool in stability theory. 

Lemma 3.1 . If f(t) and g(t) are real valued non-negative functions, 

f(t) is continuous on [tQ»“) and g(t) is integrable in every finite 

interval contained in [t ,«>) and if for some constant a > 0 

° t ~ 

(3.5) f(t) a + / f(s) g(s) ds, 

^0 

then 

t 

(3.6) f(t) exp / g(s) ds. 

t 

0 

Proof: If a > 0, then we have 



f(t) q(t) ig(t). 

a + /^ f(s)g(s)ds 
t 
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Integrating the above between the limits t^ and t yields 

t t 

In [a + / f(s)g(s)ds] - In a < / g(s)ds. 



or taking exponents, 



t t 

a + / f(s)g(s)ds < a exp / g(s)ds. 



and thus 



t 

f (t) ^ a exp / g(s)ds. 

‘o 

If a = 0, then (3.5) holds for all real numbers 6 > 0, and thus 
(3,6) holds which implies as g->-0 that f(t) = 0. This completes the 
proof. 

Utilizing Lenma 3.1, it is possible to Drove the following 
proDosition . 

00 

Proposition 3.2. If /llA(t)|| dt < then the system x' = A(t) x 

t 

0 

is stable. 

Proof: Let X(t) be the fundamental matrix of the system. Then 

from (2.6) 



so that 



X(t) = I + / A(s) X(s) ds, 
t_ 



lX(t)l| ^n + / |lA(s)|l |lX(s)|l ds. 
^0 



Aoplving Gronwall's Lemma yields 



||X(t)|| n exD / ||A(s)|| ds £ n exp / | jA(s) 

0 0 
and thus X(t) is bounded. By Theorem 3.1, the system is stable. 



ds 
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We next consider the linear homogeneous system with constant 
coefficients 

(3.7) X' = A X 

whose fundamental matrix X(t) is given by 

(3.8) X(t) = exp [(t-t^)A], 
where exn (A) is defined by 



exD (A) 



I + 



E 

n=l 



1 _ 

nl 



A 



n 



The stability of the system (3.7) is directly related to the nature 
of the roots of the characteristic polynomial of A given by 
(3.9) det (A - XI) = 0. 

Since (3.9) is an nth degree polynomial equation in x, there exists 
n roots to this equation (including repeated roots) which may be 
real or complex. Corresponding to each root is a solution of (3.7) 
and the n roots thus furnish n linearly independent solutions to 
(3.7). If X is a real root of multiplicity k, then corresponding 
to X are the k linearly independent solutions 

e*‘; te‘‘; 

while if x+iy is a complex root of multiplicity m, then so is 
x-iy and corresponding to these 2m roots are the 2m linearly indeoen- 
dent solutions 

e^^cos yt, e^^sin yt; te^^cos yt, te^^sin yt; ...; t'^'^cos yt, t’’^"^sin yt. 

Since sin t and cos t are bounded, and since for a, g >0, 

lim t e ^ = 0, we note that if the real parts of the eigenvalues of 



A (roots to (3.9)) are negative, then all solutions of (3.7) will 
converge to zero. When this is the case, we have the following 
defini tion . 
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Definition 3.4. If the real parts of the eioenvalues of A are 



negative, then we say that the characteristic polynomial of A is 
stable. 

When the characteristic polynomial of A is stable, then the 
solutions of (3,7) are of exponential order as the following lemma 
illustrates. 

Lemma 3.2 . If the characteristic nolynomial of A if stable, 
then all solutions of (3.7) satisfy 
(3.10) I |x(t)| I < Me‘“'^ t t^, 

where M and a are positive constants. 

Proof: Let A|^ = + ioj^, k = l,...,m, m_^n, be the 

distinct eigenvalues of A. Since the characteristic nolynomial of A 
is stable, there exists a > 0 such that + a < 0, k = l,...tn. 

A solution Z|^(t) corresponding to is of the form 

where bj^ is a constant vector having the entry 1 in one of its rows 
and zeros elsewhere, and where r depends on the multiplicity of A|^. 
Then since + a < 0, 

||Z|^(t)|| e“^-^0 as t-^«> 
and thus there exist constants M|^ such that 

I |Z|<(t)l I 1 for t ^ t^. 

If x(t) is any solution to (3.7), it can be expressed .as the linear 
combination 

x(t) = z.j(t). 
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where the summation is over the n basis solutions. Then 

I |x(t)| I 1 E|e.| I lz.(t)| I 

max I E e”“^ <_ M e'^^, 

where M = max 1 6^. | . 

It is now easy to show that when the characteristic oolynomial 
of A is stable, then (3.7) is an asymototical ly stable system. 

Proposition 3.3 . If the characteristic polynomial of A is stable, 
then the system x' = Ax is asymptotically stable. 

Proof: By Lemma 3.2, every solution x(t) of (3.7) satisfies 

||x(t)|| e"“^ for appropriate M and a. Since the columns of 

X(t) are solutions, where X(t) is the fundamental matrix of (3.7), 
we have 

||X(t)|| iMne'“^ t it^, 

and thus ||X(t)||->-0 as t-x» and the system is asymptotical ly stable 
by Corollary 1 to Theorem 3.1. 

Thus the stability of a linear homogeneous system with constant 
coefficients depends solely upon the nature of the eigenvalues of A. 

The question now arises as to the stability of a system (3.4) in which 
A(t) is not constant. Proposition 3.1 insures that integrability of 
A(t) is sufficient to guarantee stability in the system. However, this 
is a rather stringent requirement and is of little practical signi- 
ficance in determining the stability of a given system. It is often 
easier to determine stability in a given system by comparing it to a 
known stable system. The simolest case in which this may be done 
successfully is when A(t) is asymototic to a constant matrix B. 

The following theorem then applies 
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Theorem 3.2. If the characteristic polynomial of B is stable, 



A(t) is continuous on [0,«>), and 

oo 

/ ||A(t) - B|| dt < » . 

0 

then the system x' = A(t) x is asymptotically stable. 

Proof: The system x' = A(t) x may be written in the form 

x» = Bx + [A(t) - B] X 

and in this manner may be treated as though it were a non-homogeneous 
system with constant coefficients. The solution x(t;0,x^) may then 
be written by Proposition 2.2 and (2.13) 

t 

x(t) = X(t) X + / X(t-s) [A(s) - B] x(s) ds, 

° 0 

where X(t) is the fundamental matrix of the constant system 
x' = B X. By Lemma 3.2., there exist M,oi >0 such that 
||X(t)|| Me”“^ and thus 

llx(t)|| llxjl t ;H ||A(s)-B|l ||x(s)l| e“®ds. 

0 

Now applying Gronwalls* Lemma, 

+ t 

I |x(t) II e“^ ^ M I |Xq| I exp / M | |A(s)-B| | ds, 

0 

from which 

||x(t)|| e“^ iM ||x J| exn ||A(s)-B|| ds 

° 0 

IK, 

where K is a convenient constant. Hence ||x(t)|| £ Ke~“^ and if 
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0 3 S 



y(t;0,y^) is any other solution, 

I |x(t) - y(t)| I ^ 
and thus the system x' = A(t) x is asymptotically stable. 

Theorem 3.2, may be modified in the event A(t) does not differ 
much from B, i.e., ||A(t) - B| | is small as follows. 

Theorem 3.3 . If the characteristic oolynomial of B is stable, 
and ||A(t) - B| | ^ c on [0,oo) for sufficiently small c, then the 
system x' = A(t) x is asymptotical ly stable. 

Proof: Let x(t;0,x^) be a solution of x' = A(t) x. Then x(t) 

may be written 

t 

x(t) = X(t) X + / X(t-s)[A(s) - B] x(s) ds, 

° 0 

where X(t) is the fundamental matrix of x' = B x. Then proceeding as 
in the proof of Theorem 3.2, we have as before 

I |x(t)| I e“^ ^ M I |Xq| I exp M / | |A(s) - B| | ds, 

0 

where M and a are as previously defined in Theorem 3.2. Thus 

1 |x(t) II 1 I I^qI I (Me -a) t , 

so that if c is sufficiently small, i.e., c < then x(t)^0 as 
t^oo and since x(t) is an arbitrary solution, the system is 
asymptotically stable. 

If the characteristic polynomial of B is not stable but the solu- 
tions of x' = B X are bounded, then the following theorem holds. 

Theorem 3.4 . If all solutions of x‘ = B x are bounded on [0,«), 
A(t) is continuous on [0,oo) and 

oo 

/ I |A(t) - B| I dt < “, 

0 

then the system x' = A(t) x is stable. 
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Proof; Let X(t) be the fundamental matrix of the constant 



system x' = B x. Then there exists M > 0 such that | |X(t)| | <_ M 

for t ^ 0. If x(t;0,x^) is a solution of x' = A(t) x, then 

t 

x(t) = X(t) X + / X(t-s)[A(s) - B] x(s) ds. 
° 0 



Thus 



l|x(t).|| < M ||x II + ||A(s) - B|| ||x(s)|| ds, 

° 0 



and applying Gronwall's Lemma, 

t 

||x(t)|| ||x II exp f M ||A(s) - B| | ds, 

° 0 



from which 



||x(t)|| M l|x II exp f M ||A(s) - B|| ds 
° 0 

^ K < oo 

where K is a suitable constant. The above shows that all solutions are 
bounded and hence the system x' = A(t) x is stable by Theorem 3.1. 

Theorems 3.2, 3.3, and 3.4 insure that if ||A(t) - B|| is well 
behaved, then the solutions of x' = A(t) x will resemble those of the 
constant system x' = B x as t^ when A(t) is asymptotic to B. In 
view of this, we might expect these results to generalize to the case 
when B is not constant. This is not always the case though it is 
possible to effect a comparison between two non-constant systems under 
more restrictive hypotheses. Toward this end, we prove the following 
1 emma . 
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and 



Lemma 3.3 . If the system x' = A(t) x is stable on 



t 

lim inf / trA(s) ds > - « , 

t-><» t 

0 

then X~^(t) is bounded on [t^,°°), where X(t) is the fundamental matrix 
of the system. 

Proof: Let X(t) be the fundamental matrix of the system. 



Then 



x-’(t) 



adj X(t) 
det X(t) 



and since x' = A(t) x is stable, it follows that X(t) is bounded for 

t > t and this implies that there exists L > 0 such that 
— 0 

1 1 ad j X ( t ) 1 1 ^ L for t ^ t^ . 

Then this relation and (2.8) yield 

||X-’(t)|l < ^ 

exp / tr A(s) ds 



and since 



t 

lim inf / tr A(s) ds > - «>, 

t-^oo t 

0 

it follows that 



exp / tr A(s) ds ^ k > 0, 
^0 

so that 

l|X“^(t)|| for t ^t^ 
and thus X”^(t) is bounded. 
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We can now state the fol lowing comoarison theorem for two non- 
constant systems. 

Theorem 3.5 . If the system x' = B(t) x is stable on [t^,«>). 



t 

lim inf / tr B(s) ds > - «, 
t->«> t 

0 

A(t) is continuous on and 

oo 

/ I |A(t) - B(t)| I dt < 00, 

‘o 

then the system x* = A(t) x is stable. 

Proof: Let X(t) be the fundamental matrix of the system 

X* = B(t) X. Then the hypotheses imply that there exists M > 0 
such that for t ^ t^, ] lX(t) | | and | |X“^ (t) | | ^M. 

Then if x(t;t^,x^) is any solution to x' = A(t) x, x(t) may be written 



x(t) = X(t) X + / X(t) X"^(s) [A(s) - B(s)]x(s)ds, 



and utilizing the above we obtain 



||x(t)|| ||x II * ! l|A(s) - B(s)|| ||x(s)||ds. 

‘o 

Applying Gronwall's Lemma yields 



so that 



|x(t)| I ^ M I |x^| I exp 




1 |A(s) - B(s)| 1 ds. 



|lx(t)|| l|x II exD / ||A{s) - B(s)|lds 

1 K, 

where K is a convenient constant. Thus all solutions of x* = A(t) x 
are bounded and hence stable by Theorem 3.1. 
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Theorem 3.5 gives one set of conditions under which a comoarison 
may be effected between two non-constant linear homogeneous systems 
and is a well known result. Another set of conditions which have 
been formulated by the author and which are also sufficient to effect 
a comparison are stated below in the following theorem. 

Theorem 3.6 . If the system x' = B(t) x is stable on [t^,®), 

A(t) is continuous on [t^,®) and 



where X(t) is the fundamental matrix of x' = B(t) x, then the 
system x' = A(t) x is stable on [t^,«>). 

Proof: Let X(t) be the fundamental matrix of x' = B(t) x. Then 

if x(t;t^,x^) is any solution to x' = A(t) x, x(t) may be written 



f I IX"'' (t) [A(t) - B(t)] X(t)| I dt < « 
t 



t -| 

x(t) = X(t) X + / X(t) X (s) [A(s)-B(s)]x(s) ds. 
° t 



0 



From this we obtain 



X''(t)x(t) 




^ -1 

/ X (s)[A(s) - B(s)]x(s) ds 
t_ 




x(s)ds 



and thus 



||X-'(t)x(t)|| i llxjl + 



t 
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so that applying Gronwall's Lemma yields 



l|X‘''(t)x(t)|| il|x„|| exp / ||X-'(s)[A(s)-B(s)] 

‘o 

X(s)|| ds. 

Therefore, 

I |X"’(t)x(t)| I ^ I |Xj|| I exp / I |X”’(s)[A(s)-B(s)] 

° X(s)|| ds 

^ K < “ 

which shows that X~^(t)x(t) is bounded. 

However, 

l|x(t)|| = II X(t)X-'(t)x(t)|| i ||X(t)|| ||X-'(t) 

x(t)|| 

SO that ||x(t)|| K ||X(t)|| and since x' = B(t) x is stable, X(t) is 
bounded and thus x(t) is bounded. Hence all solutions to x' = A(t) x 
are bounded and thus stable by Theorem 3.1. 

Theorems 3.2 through 3.5 provide conditions under which the 
stability of a given system may be determined by a comnarison with 
a known stable system. Other conditions under which similar 
theorems can be established have been given by Bellman [1], and 
Levinson [6], as well as numerous other authors. 
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IV. NON-HOMOGENEOUS LINEAR SYSTEMS 



The systems considered here are of the form 
(4.1) x' = A(t) X + b(t ) 

where A(t) and b(t) are continuous on [t^,“). 

Recall that for the homogeneous system, b(t) e 0, boundedness 
and stability of the system are equivalent. This is no longer true 
for the system (4.1) as the following example illustrates. 

Example 4.1 . The one-dimensional system 
x' = 2t on [0,oo) 
has solutions x(t'»0,x^) of the form 

x(t) = + t^ 

which are stable yet unbounded. 

The above example shows that it is possible for a non-homogeneous 

to be stable but unbounded. The reverse behavior, however, is not 

possible as the following theorem shows. 

Theorem 4.1 . If the system x' = A(t) x + b(t), A(t) and b(t) 

continuous on [t ,«>), is bounded on [t ,») then it is stable on 
0 0 

Proof: Let X(t) be the fundamental matrix of the homogeneous 
system x' = A(t) x. If (4.1) is bounded, then in particular, the 
solution x(t;t^,0) given by 

x(t) = X(t) / X '(s) b(s) ds 
t 

0 

is bounded. Thus there exists M.j> 0 such that ||x(t)|| M^ ^ 
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If y(t',t^,yQ) is any other solution to (4.1), then 

t 1 

y(t) e X(t)y + f X(t)x (s) b(s) ds 
^0 

and there exists such that 

||y(t)l| = I |X(t) y^ + x(t)| I ^ M 2 for t ^ t^ 

Then 

l|X(t) yjl - ||x(t)|| 1 M 2 , 

so that 

l|x(t) yjl iM t ^t^, 

where M = M, + M^. Therefore if w(t;t ,w ) and v(t;t ,v ) are any 
\ c 0 0 0 0 

two solutions of (4.1) and if e > 0, then for t > t 

— 0 

I |w(t) - v(t)| I = ||x(t) (Wq-Vq)II < M ||Wq-Vq 

so that if ||Wq-Vq|| < then | |w(t)-v(t) | | <e ^ 
which shows that (4.1) is stable. 

Theorem 4.1 admits the following partial converse. 

Theorem 4.2 . If the system (4.1) is stable and one solution is 
bounded, then all solutions are bounded. 

Proof: Let X(t) be the fundamental matrix of the homogeneous 
system x' = A(t) x. Once again the first part of the proof consists 
of showing that X(t) is bounded. 

If (4.1) is stable, then in particular, the solution x(t;t^,0) 
given by 

f -1 

x(t) = / X(t)X (s)b(s) ds 

^0 

is stable. Therefore, if e > 0, there exists 6 > 0 such that 

for t > t , 

— 0 
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1 |w(t)-x(t) 11 = 1 I X(t)w^| I < e 

whenever w(t;t ,w ) is a solution with llw || < 6. Letting 
Wq = where e^ is the vector whose jth component is el = 6^^, we 

obtain from the above 

ll♦i(t)ll ' |£ 

where (fi^(t) is the ith column of X(t). Then since this is valid for 
i = 1 ,...,n, 

I |X(t)| I < ^ = M. 

Now let z(t;t^,z^) be a bounded solution and y(t;t^,y^) any arbitrary 
solution of (4.1). Then since y(t) = z(t) + [y(t) - z(t)], 

I ly(t)| I ^ I |z(t)| I + I |X(t)[y^-z^]| I 

ll |z(t)| I + M I ly^ - z^l I 

which shows that y(t) is bounded. Hence the system (4.1) is bounded. 

Theorems 4.1 and 4.2 show that in a non-homogeneous linear system, 
boundedness is a stronger property than stability. For this reason, 
it is usually easier to determine stability in a system (4.1) than it 
is boundedness. In fact, the following proposition shows that the 
stability of (4.1) may be determined by simnly examing its homo- 
geneous part. 

Proposition 4.1 . If the system x' = A(t) x is (asymptotically) 
stable on [t^,“), then so is the system x' = A(t) x + b(t), where 
b(t) is continuous on [t^,~). 

Proof: If x(t;t^,x^) and y(t;t^,yp) are two solutions of 

x' = A(t) X + b(t), then the difference [x(t)-y(t)] satisfies the 
homogeneous system x' = A(t) x. Thus since the trivial solution 
z(t) = 0 of the homogeneous system is stable, given e > 0, there 
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exists 6 > 0 such that if | | (x^ ■ Yq) ■ 0|| <5 then 

I |x(t) - y(t) - 0| I < e t ^ t^, 

which is precisely the statement that the non-homogeneous system is 
stable. The proof is the same if the homogeneous system is, in 
addition, asymptotically stable. 

If A(t) is a constant matrix A, then (4.1) becomes 
(4.2) x' = A X + b(t) 

and boundedness is easier to determine. If the characteristic 
polynomial of A is stable, then the following theorem holds. 

Theorem 4.3 . If the characteristic polynomial of A is stable, 
b(t) is continuous and ||b(t)|| ^ K, for t ^ 0, then the system (4.2) 
is bounded. 

Proof: Let X(t) be the fundamental matrix of the constant 

system x' = A x. By Lemma 3.2, there exist M and a such that 
I |X(t) 1 1 Me'“^ for t ^ 0. If x(t;t^,x^) is a solution of (4.2) 

then from 



t 

x(t) = X(t)x + f X(t-s)b(s) ds, 
° 0 



we deduce successively 



l|x(t)|| iM llx^ll + 



/ ds 

0 



- M 



. ^ (l-e-“t) 

a 



< M ||x J| + ^ 

which shows that x(t) is bounded. Thus the system (4.2) is bounded, 
A similar theorem holds when the solutions of x' = A x are 
bounded, though the forcing term, b(t), must be integrable. 
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Theorem 4.4. If the constant system x' = A x is bounded on 



[0,«), b(t) is continuous on [0,“) and 

oo 

/ I ] b(t) I I dt < «> , 

0 

then the system x' = A x + b(t) is bounded. 

Proof: Let X(t) be the fundamental matrix of x' = A x. Then 

there exists M > 0 such that ||X(t)|| If x(t;0,x^) is 

a solution to x' = A x + b(t), then from 

t 

x(t) = X(t) X + / X(t-s)b(s) ds, 

° 0 

we deduce successively 

t 

||x(t)|| < M I |x II + / M ||b(s)|| ds 
° 0 

oo 

^M||x|| + / M|| b(s) I I ds 

° 0 

^ K < “ 

where K is a convenient constant. Thus the system x' = A x + b(t) 
is bounded. 

The boundedness of a non-homoqeneous system x' = B(t) x + b(t) 
when B(t) is asymptotic to a constant matrix A may be determined by 
effecting a comoarison between it and the constant svstem x' = A x. 
The following theorem gives one instance when this may be done 
successful ly . 

Theorem 4.5 . If the characteristic polynomial of A is stable, 

b(t) is continuous and bounded on [0,®), B(t) is continuous on [0,°°) 

and /°°||B(t) - A| | dt < “, then the system x' = B(t) x + b(t) is 
0 

bounded and asymptotically stable. 
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Proof: Let X(t) be the fundamental matrix of the constant system 

x‘ = A X. Then there exist M, a > 0 such that ||X(t)|| Me”“^ 

for t ^ 0. If x(t;0,x^) is a solution to x' = B(t) x + b(t), 
then as before, 

t 

x(t) = X(t)x + ; X(t-s)[(B(s)-A)x(s) + b(s)] ds, 

° 0 



from which 



l|x(t)|| i l|X(t)|| ||x II + ||X(t-s)|| 

° 0 

(||B(s)-A|| ||x(s)|| + ||b{s)||) ds. 

Then since | | b(t) | | ^ N 

||x(t)|| e“* < M ||x J| + / Me-“^||R(s)-A|| ||x(s)||ds 
t 

+ f MNe'“^ ds 
0 

< M ||x l| + / Me-“"||x(s)|| ||B(s)-A||ds 
° 0 

00 

+ f MNe"“^ ds 
0 

and thus 

||x(t)|| e“\ L+ / Me-“®||x(s)|| ||B(s)-A|| ds, 

0 

where L is a convenient constant. Aonlyinq Gronwall's Lemma yields 

I |x(t) I |e“^ £ L exD f M | |B(s)-A| | ds, 

0 

so that 



||x(t)l|e“^ < L exD f M | |B(s)-A| | ds < K < «> 

0 
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where K is another constant. Thus for t ^ 0, |•|x{t)|| ^ Ke”^^ ^ K, 
which shows that x(t) and hence the system x' = B(t) x + b(t) is 
bounded. 

In addition, if y(t;0,y^) is any other solution to this system, 

I |x(t) - y(t) I I <_ I |x(t) I I + I |y(t) I I ^ 2Ke"“^ 
which shows that | |x(t) - y(t) | l->0 as t->~ and the system 
x' = B(t) X + b(t) is also asymptotically stable. 

If the characteristic polynomial of A is not assumed stable, 
but the system x' = A x is bounded, then a similar theorem holds. 

Theorem 4.6 . If the system x* = A x is bounded on [O,"), b{t) 
and B(t) are continuous on [0,“), 

oo 

/ I I b(s) II ds < ■», and 

0 

oo 

/ 1 I B(s) - A| I ds < “ , 

0 

then the system x' = B(t)x + b(t) is bounded. 

Proof; Let X(t) be the fundamental matrix of the constant 
system x' = A x. Then there exists M > 0 such that j |X(t)| | M for 
t ^ 0. If x(t;0,x^) is any solution to x‘ = B(t) x + b(t), then from 

t 

x(t) = X(t)x + / X(t-s)[(B(s)-A)x{s)+b(s)] ds, 

° 0 



we have successively 



1 M I lx I I + / M 1 |B(s)-Al I I |x(s)l 1 ds 
° 0 



+ / M I I b(s) I I ds 

0 



< K + / M| lB(s)-A| I I |x(s)| I ds, 

0 

where K is a suitable constant. Applying Gronwall's Lemma yields 



|x(t)| I < K exp / M I |B(s)-A| | ds 
0 



< K exp / M I |B(s)-A| | ds < L < «> 

0 

where L is another constant. Thus x(t) and hence the system 
x' = B(t) X + b(t) are bounded. 

In general, if A(t) is not constant, Theorems 4.3 through 4.6 
need no longer apply. Thus if a homogeneous system x' = A(t) x 
is bounded (stable) and 

/ I |b(t) II dt < -, 

0 

it need not follow that the non-homogeneous system x' = A(t) x + b(t) 
is bounded. The same applies if the homogeneous system is asymptotically 
stable. A similar problem was encountered in attempting to effect a 
stability comnarison between two non-constant linear homogeneous systems. 
There it was noted that additional hypotheses were required to insure 
the stability of the system x' = B(t) x when B(t) is asymototic to a 
non-constant matrix A(t) and the system x' = A(t) x is stable. One 
reason for this is that stability, as defined in the sense of Lyaounov, 
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is quite delicate and often may not be maintained even when small 
oerturbations are introduced. This leads to the concepts of uniform 
and restrictive stability, which are discussed in section V. 
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V. UNIFORM AND RESTRICTIVE STABILITY 



Since Lyapunov stability does not suffice to insure certain 
boundedness and comparison properties, we are led to a consideration 
of stronger forms of stability. One of these is uniform stability, 
defined as follows. 

Definition 5.1 . Let x(t;t^,x^) be a solution of the system 
x' = f(t,x), where f(t,x) is defined and continuous in the region 
F = {(t,x): 1 t < 00 , I |x| I < 5 }. Then x(t) is said to be uniformly 

stable if it satisfies (3.2i,ii, and iii) and in addition (5.1) 
given e > 0 there exists 6 = 6(e,f) > 0 such that every solution 
y(t;t^,yQ) with ||x(t^) -y(t.j)|| < $ for some t.j ^ t^ satisfies 

I |x(t) - y(t) II < e for t ^ t^ . 

Comnaring the definitions, we note that uniform stability is 
similar to Lyapunov stability. However, there is an important 
distinction between the two. A solution x(t) is stable in the sense 
of Lyapunov if any other solution which is initially close to x(t) 
remains in an e-neighborhood of x(t) for t ^ t^; at later times, 
i.e., t > t^ ; other solutions may come arbitrarily close but need 
not remain there. If x(t) is uniformly stable, once a solution enters 
a prescribed small neighborhood of x(t) for any t-j ^t^> it must remain 
there for all t ^ t-| . Thus, in uniform stability, 6 no longer 
depends on t^. 

The above remarks show that a uniformly stable solution is stable. 

If all solutions to a system are uniformly stable, we say that the 
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system is uniformly stable. For linear homogeneous systems of the 
form 

(5.2) x' = A(t) X, 

if one solution to the system is uniformly stable, then the system 
is uniformly stable. Uniform stability for a system (5.2) may be 
characterized by the behavior of the fundamental matrix of the 
system as the following theorem shows. 

Theorem 5.1 . The system x' = A(t) x is uniformly stable iff 
there exists M > 0 such that | |X(t)X"^ (s) | | ^ M t ^ s ^ t^, where 
X(t) is the fundamental matrix of the system. 

Proof: Assume that the system x' = A(t) x is uniformly stable, 

so that the trivial solution z(t) e 0 is uniformly stable. Then 
if e > 0 there exists 6 > 0 such that if ||x(t^)|| < & for some 
t-| ^ tpj then I |x(t) 1 1 < e for t ^ t.| . Let y(t;t^,y^) be a solution 
of x‘ = A(t) X. Then y(t) may be written 

y(t) = X(t) X"''(t^) y(t^). 

By choosing y(t) to be that solution satisfying y(t,) = \ where 

i 5 

e is as in the proof of Theorem 3.1, we have ||y(ti)|| = ^< 6 

so that 

|| |X(t)X“''(t^)e‘'| I < £ t^t^, 

or eguivalently, 

||z'(t,t^)||< t ^t^ 

where Z^(t,t-|) is the ith column of the matrix X(t)X"^ (t-j ) . Repeating 
the above orocedure for i = l,...,n yields 

||X(t)X'''(t^)|| < ^ t ^t^. 
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Finally, since this procedure is valid for any t-j ^ and the 

bound is indeoendent of t-| , we obtain 

I |X(t)X‘"'(t^}| I < ^ t > t^ ^ t^. 

Conversely, if j |X(t) X"^(t-|)| | ^ M for t ^ t-| ^ t^, and 

x(t;t^,XQ) and y(t>tQ»yQ) any two solutions of x' = A(t) x, 

then for t > t, > t 
— 1 — 0 

I |x(t)-y(t)| I = I |X(t)X"^(t^)[xCt^)-y(t^)]| I 
1 M I |x(t^} - y(t^)| I 

so that if e > 0, I |x(t^ ) - y(t-| ) 1 1 < ^ 1 

||x(t) - y(t)|| < e t ^ t^ which shows that the system is uniformly 
stable. 

From Theorem 5.1, we obtain the following simple proposition. 
Proposition 5.1 . If the system x' = A x, where A is constant, 
is stable on [0,oo), then it is uniformly stable on [0,oo). 

Proof: Let X(t) be the fundamental matrix of x' = A x. Then 
there exists M > 0 such that ||X(t)|| for t ^ 0. Then by (2.13), 

||X(t-s)|l = I |X(t) X'^(s)| I ^ M for t ^ s ^ 0 

and the system is thus uniformly stable. 

Recall that in sections III and IV, several boundedness and 
comparison properties which hold when A(t) = A is constant, need not 
apply when A(t) is not constant. However, Proposition 5.1 shows 
that in the case A(t) is constant, the systems considered are actually 
uniformly stable, rather than just stable. This leads to the question 
of whether or not uniform stability of the non-constant system 
x' = A(t) X is sufficient to insure these boundedness and comparison 
properties. The following theorems provide an affirmative answer 
to this question. 
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Theorem 5.2. If the system x‘ = A(t) x is uniformly stable on 



[t ,«>), B(t) is continuous on [t and ^./°°| | B(t)-A(t) | | dt 
0 ° ^0 

then the system x' = B(t) x is uniformly stable. 

Proof: Let X(t) be the fundamental matrix of the system 

x' = A(t) X. Then ||X(t) X"^(s)|| for t >_ s ^ t^. Let e > 0 

and fix s ^ t^. If x(t;t^,x^) and y(t;t^,y^) are two solutions of 
the system x‘ = B(t)x, then since B(t) = A(t) + [B(t) - A(t)], 
we have for t >_ s, 

x(t)-y(t) = X(t)X"^ (s)[x(s)-y(s)] 

+ / X(t)X ^ (T)[B(T)-A(T)][x(x)-y(T)] dx 
s 

and thus 

l|x(t)-y(t)|| < ||X(t)X-'(s)|| ||x{s)-y(s)|| 

+ / II X(t)X ^(t)|| I I B(t)-A(t) I I I |x(t )- 

y(x) I |dx 



^ M I |x(s)-y(s) I I 



+ / M I I B(x)-A(x ) I I I |xCx)-y(x) 1 1 dx. 

t 

0 

Applying Gronwall's Lemma yields 



from which 



x(t)-y(t)|| I |x(s)-y(s)| I exn 



/ M ||B(x)-A(x)|| 
t 



00 

I |x(t)-y(t)| I I |x(s)-y(s)| I exp / M ||B(x)- 

^0 

A(x) I I dx . 
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Therefore , 



I |x(t)-y(t)| I 1 L I |x(s)-y(s)| | 

where L is a convenient constant. The above relation is valid 
for t ^ s ^ t^, so that if | |x(s)-y(s) || < ^ ^ L 

I |x(t)-y(t) I I < e for t ^s which shows that the system x' = B(t) x 
is uniformly stable. 

Corollary 1 . Under the same hypotheses as Theorem 5.2, the 
system x' = B(t) x + b(t), where b(t) is continuous on [tQ,“) is 
uniformly stable on [t^,<»). 

Proof: The proof follows from the fact that under the hypotheses, 

the system x' = B(t) x is uniformly stable and that if x(t) and 
y(t) are two solutions to the non-hompgeneous system x* = B(t) x + b(t), 
then the difference x(t) - y(t) is a solution to the non-homogeneous 
system x' = B(t) x. 

When the homogeneous system is uniformly stable, then 

boundedness of the non-homogeneous system may be determined as 

for the constant case considered earlier in Section IV. 

Theorem 5.3 . If the system x' = A(t) x is uniformly stable 

on [t ,<”), b(t) is continuous on [t ,«) and /'f|b(s)|| ds < ”, 

0 ° 0 

then the system x' = A(t) x + b(t) is bounded on [t^,“). 

Proof; Let x(t;t^,x^) be a solution to x* = A(t) x + b(t). 

Then 

t 

x(t) = X(t)x + / X(t)X (s)b(s) ds, 

^0 

where X(t) is the fundamental matrix of the homogeneous system 
x' = A(t) X. Since llX(t) X~^(s)|| for t L s ^ t^, 
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and thus 



t 

x(t)|l ||x || + / H ||b(s)|| ds 
° t 

0 

oo 

|x(t)|| ||x II + / M ||b(s)|| ds ^L< « 

where L is a convenient constant. Thus the system x' = A(t) x + b(t) 
is bounded. 

It is also possible to determine boundedness of a non-homoqeneous 
system through comparison with a uniformly stable system. 

Proposition 5.2 . If the system x‘ = A(t) x is uniformly stable 

CX3 

‘on [t , B(t) and b(t) are continuous on [t ,»), /||t»(t)|| dt < «>, 

0 0 ^ 

0 

oo 

and /||B(t) - A(t)|| dt < then the system x* = B(t) x + b(t) 
t 

0 

is bounded on [t ,“). 

0 

Proof: By Theorem 5.2, the hypotheses imply that the homo- 

geneous system x' = B(t) x is uniformly stable on [t^joo). The 
result then follows from Theorem 5.3. 

Thus certain boundedness and comparison properties established 
for constant systems may be generalized to non-constant systems 
which are uniformly stable. 

We now introduce the concept of restrictive stability. 

Definition 5.2 . The linear homogeneous system x' = A(t) x 
is said to be restrictively stable if it is stable and its adjoint 
system x' = -A(t)^ x is stable. 

We note that restrictive stability applies to a system, rather 
than a particular solution and is defined for linear homogeneous 
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systems only. The following theorem characterizes restrictive 
stability in terms of the behavior of the fundamental matrix of the 
system. 

Theorem 5.4 . The system x' = A(t) x is restrictively stable 
on iff there exists M > 0 such that ||X(t)|| and 

1 1 X"^ (t) 1 1_< M for t ^ where X(t) is the fundamental matrix 
of the system x' = A(t)x. 

Proof: The result follows by Theorem 3.1, noting that [X~^(t)]^ 

is the fundamental matrix of the adjoint system x' = -A(t)^x. 

From Theorem 5.4, we see that a restrictively stable system 

x' = A(t) X is uniformly stable since 

||X(t)X-'(s)|| < Ilx(t)|l ||X-'(s)|| < 

T T 

but not necessarily conversely. Also, since -[-A(t) ] = A(t), 

it follows that if x' = A(t) x is restrictively stable, then its 
adjoint system is restrictively stable. The following oropositions 
give conditions which suffice for a system to be restrictively 
stable. 

c» 

Proposition 5.3 . If / ||A(t)|| dt < «, and A(t) is continuous 

‘o 

on [t^,«), then the system x' = A(t) x is restrictively stable on 

[tp,”) . 

Proof; Proposition 3.2 shows that under the hypotheses, the 
system x' = A(t) x is stable. Since ||A(t)|| = 1| -A(t)^||, it 

follows that the adjoint system is stable. 

Proposition 5.4 . If the system x' = A(t) x is stable on [t^,") 
and A(t) is continuous on [t ,■»), and lim inf tr A(s) ds > 

° t->=o t 

0 

then it is restrictively stable. 
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Proof: Let X(t) be the fundamental matrix of x‘ = A(t) x. Then 

X(t) is bounded for t ^ t^ since the system is stable. By Lemma 3.3, 

lim inf } tr A(s) ds > implies that X”^(t) is bounded for 
t-x» t 

0 

t ^ t^. Hence the system x' = A(t) x is restrictively stable by 
Theorem 5.4. 

Proposition 5.5 . If the adjoint system x* = -A(t)^x is stable 
on [t ,«>), A(t) is continuous on [t ,«>), and lim suo tr A(s) ds < 

° t->co t„ 

then the system x' = A(t) x is restrictively stable on [t^,«). 

Proof; Let X(t) be the fundamental matrix of x' = A(t) x. Then 
since the adjoint system is stable, X“^(t) is bounded for t ^ t^. 

From the relation X(t) = det X(t) adj X”^(t), we obtain by (2.8) 

t 

X(t) = exp / tr A(s) ds adj X (t) 
t 

0 

and since X"^(t) is bounded, adj X”^(t) is bounded and thus 

t 

I |X(t)| I ^ M exp r tr ACs) ds 
^0 

t 

lim sup / tr A(s) ds < L < <» 
t->«> t 

0 

where M and L are suitable constants. Hence X(t) is bounded and 
the system x' = A(t) x is restrictively stable by Theorem 5.4. 

An interesting result concerning the well known second order 
Hill's Equation 

(5.3) u" + a(t) u = 0 on [O,-), 

where a(t) is real valued and continuous on [0,“) and u(t) is 

a scalar valued function, may be obtained through the use of 
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stability theory. The vector analogue of (5.3) is the two-dimensional 
system 



(5.4) y' = z 

z' = -a(t) y 



which is in the form (5.2) with x = and A(t) 



0 



-a(t) 0^ 



We note that, since tr A(t) = 0, if the system (5.4) is stable, 
then it will be restrictively stable, and thus uniformly stable. 



It is now possible to give the following result: If /la(t)| dt < «>, 



then there exist unbounded solutions of (5.3). This may be 
proved using stability theory as follows. Assume that all solutions 
to (5.3) are bounded. Then (5.4) is a stable system and hence 
uniformly stable by the above remarks. Theorem 5.2 then asserts 
that if 




then the system x' = B(t) x is uniformly stable and hence bounded 

OO 00 

since / |lB(t) - A(t)l| dt = / |a(t)| dt < «>. But this is a 

0 0 

contradiction since 



is a solution to x' = B(t)x. Thus all solutions to (5.3) are bounded. 

Restrictive stability may also be characterized in terms of 
the reducibility of the system. Toward this end, we have the 
following definition. 



0 



x(t) 
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Definition 5,3. Tne system x' = A(t) x is said to be reducible 



if there exists a matrix L(t) whose elements are differentiable 
functions which is bounded together with L”^(t) on [t^,") such 
that 

(5.5) L'^(t) A(t) L(t) - L"^(t) L'(t) = B 

where B is a constant matrix. If B = 0, the zero matrix, then 

we say that the system is reducible to zero. 

The relation (5.5) is often exnressed by saying that A(t) is 
kinematically similar to the constant matrix B. If both L(t) and 
A(t) are constant, then (5.5) becomes 

(5.6) L"^ A L = B 

which is the ordinary similarity relation. 

If a linear homogeneous system is reducible, then the trans- 
formation x = L(t) y reduces the system x* = A(t) x to the constant 
system y' = B y, while if x' = A(t) x is reducible to zero, then 
this transformation yields the system y' =0. 

The following theorem characterizes restrictive stability in 
terms of the reducibility of the system. 

Theorem 5.5 . The system x' = A(t) x is restrictively stable on 
[Tq,") iff it is reducible to zero on [t^,oo). 

Proof: If x' = A(t) x is restrictively stable, then there 

exists M > 0 such that | |X(t) M ^ M and MX"^ (t) | | ^ M for t ^ t^, 
where X(t) is the fundamental matrix of the system. Then since 
X(t) satisfies X(t) = A(t) X, it follows that 

X"^(t) A(t) X(t) - X‘^(t) X'(t) = 0 
and that the system x' = A(t) x is reducible to zero by taking 
L(t) = X(t). 
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Conversely, assume that (5.5) is satisfied for B = 0. Then 
from (5.5), L(t) satisfies 

L' = A(t) L 

and since L(t) is non-singular, it follows that L(t) is a basis 
matrix for the system x' = A(t) x. Hence X(t) = L(t) W, where X(t) 
is the fundamental matrix and W is a constant non-singular matrix. 

Thus ||X(t)|l i ||L(t)l| !|W|| and ||X-'(t)|| ^ ||L-’(t)|| HW’ll which 
shows that both X(t) and X" (t) are bounded. Thus the system 
x' = A(t) X is restrictively stable by Theorem 5.4. 

Uniform stability is also related to reducibility as the 
following theorem illustrates. 

Theorem 5.6 . If the system x' = A(t) x is stable and reducible 
on [0,“), then it is uniformly stable on [0,«>). 

Proof: Let X(t) be the fundamental matrix of x* = A(t) x and 

let L(t) be as in Definition 5.3. Then the transformation 
y(t) = L"^(t) x(t) reduces x' = A(t) x to the system y^ = B y, 
where B is constant. Also, since L"^(t) X(t) is non-singular and 
satisfies 

(L"'' X)' = B (L'”' X), 

it follows that L"^(t) X(t) is a basis matrix for y' = B y. Then 
Y(t) = L"^(t) X(t) W, where Y(t) is the fundamental matrix of 
y' = B y and W is a constant non-singular matrix. Since x' = A(t) x 
is stable, X(t) is bounded and by definition, both L(t) and L'^t) 
are bounded. Thus there exists M > 0 such that [|X(t)|| ^ M, 

||L(t)|| ^ M, and ||L~^(t)|| M, so thatl|Y(t)|| l|WM» which 

shows that Y(t) is bounded and thus that the system y' = B y is stable. 
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Then by Proposition 5.1, y' = B y is uniformly stable and thus there 
exists N > 0 such that | |Y(t) Y"^ (s) | | N for t ^ s ^ 0. Then 
s i nee 

X(t) = L(t) Y(t) W"^ 
we have for t >_ s ^ 0 

l|X(t)X-'(s)|| < ||L(t)|| ||Y(t)Y-'(s)|| ||L-'(t)|| 
^ N 

and thus x' = A(t) x is uniformly stable by Theorem 5.1. 

The above remarks relate the connection between uniform and 
restrictive stability in a system and the reducibility of that 
system. Also, we have seen that uniformly stable systems and 
restrictively stable systems satisfy certain boundedness and 
comparison properties that need not necessarily hold for stable 
systems. In Section VI, we shall show that a similar behavior 
is exhibited when the systems considered are non-linear. 
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VI. NONLINEAR SYSTEMS 



The systems considered here are non-linear mixed systems of 
the form 

(6.1 ) x' = A(t) X + f (t,x) 

where A(t) is continuous on [0,co) and f(t,x) satisfies certain 
aonrooriate conditions. Such systems are extremely important 
since the general first order normal system 

(6.2) x' = g(t,x) 

can often be put in the form (6.1) by taking f(t,x) = g(t,x) - A(t) 
where A(t) depends on the form of g(t,x). Conditions on f(t,x) 
are thus sought which insure stability in the non-linear system 
(6.1) when it is known that the linear portion 

(6.3) x‘ = A(t) X 

is stable. When A(t) = A is constant, then (6.1) becomes 

(6.4) x' = A X + f(t,x) 

and stability is easier to determine. 

Consider the following conditions on f(t,x). 

(6.5i) f(t,x) is continuous in the region 

r = {(t,x): 0 t < ~, I |x| 1 < o 

(6.5ii) [|f(t,x)|| = o(||x||) uniformly in t as 

||x||->0, i.e., there exists a function u(r) > 0 
such that u(r)->0 as r-^-0 and 
||f(t,x)l| 1 ||x|l u(||x||). 
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Condition (6-5i) assures existence of solutions for small 
||x|| but not uniqueness. Note also that (6.5ii) implies that 
z(t) E 0 is a solution to (6.4). 

When f(t,x) satisfies (6.5i,ii), we might heuristical ly 
reason that since f(t,x) is small when ||x|| is near zero, the 
system (6.4) might be expected to behave in a manner similar 
to x' = A X. Such reasoning is correct as the following theorem 
illustrates. 

Theorem 6.1 . If the characteristic oolynomial of A is stable 
and f(t,x) satisfies (6.5i, ii) then the solution z(t) e 0 of 
the system x' = A x + f(t,x) is asymptotically stable. 

Proof: We first show that solutions x(t;0,x^) exist and are 

defined on [0,“) when I |x^| | is small. Let X(t) be the fundamental 
matrix of the constant system x' = A x. Then X(t) satisfies for t ^0 
||X(t)|| M e'“^ for suitable M and a. The solution x(t;0,x^) of 
x' = A X + f(t,x) satisfies (where it exists) 

t 

x(t) = X(t) x^ + / X(t-s) f(s,x(s)) ds, 

° 0 



so that 



t t ^ 

||x(t)||e“^ llx^ll + f Me“^||f(s,x(s))|| ds. 

0 

The above inequality is valid for t [0,t-|) where |lx(t)|| < ^for 
t [0,t^), assuming Mx^ll <5- Then from (6.5ii), given e > 0, there 
exists 6 > 0 such that for t ^ 0, 

I |f(t,x)| I ^ el |x| I 

whenever | |x| | < 6. 
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If we take ||XqII < <5, then by continuity of x(t), there exists 
t 2 > 0 such that | |x(t)| | < 6 for 0 t < t 2 - Therefore for 
0 ^ t < T, where T = min (t-|,t 2 ), we have 

||x(t)|| e“^ M I |x I I + Me / e“^||x(s)|| ds. 

° 0 

Applying Gronwall's Lemma to the above yields 

• t t 

I |x(t) I |e“^ < M| |x II exp / Me ds 
° 0 



and thus for 0 t < T 

I |x(t)| I 1 M I |Xq| I exp(Me-a)t. 
Choosing ^ yields 

I |x(t)| I ^ M I |x^| I 0 £ t ^ T 
so that if I |x^| ! < 2 ^ . 



I |x(t) ||<|-<6 0_^t<_T. 

Since f(t,x) is defined and continuous for t ^0 and ||x|| < C. 
the above relation imolies that the solution x(t;0,x^) can be 
extended interval by interval oreservinq the above bound. Thus, 
any solution x(ti0,x^) with I 1x^1 | < exists on [0,~) and satisfies 

||x(t)|| < 2 " for t ^ 0. Since 6 can be taken arbitrarily small, 
this implies that the solution z(t) h 0 of x' = A x + f(t,x) 
is stable. Furthermore, 

I |x(t) I I 1 M I |x^| I exD (Me-a)t 

where e < ^ implies that ||x(t)||^0 as t^“ and thus the trivial 
solution is asymptotically stable. 



53 



Coddington and Levinson [3] show that the hypotheses of Theorem 6.1 
may be weakened somewhat. Instead of requiring that f(t,x) satisfy 
(6.5ii) it is sufficient to assume instead that for some K > 0, 

I |f(t,x)| I ^ K I |x| I t ^ 0 

for all small ||x||, and that given e > 0 there exist 6,T such that 

||f(t.x)|| < s ||X|| 
for I I X I I < 6 and t ^ T. 

Theorem 6.1 may then be applied to a system of the form 
(6.6) x' = A X + B(t) X + g(t,x) 

where B(t) ->0 as t->-“ and ||g(t,x)|| = o(||x||). 

Other conditions on f(t,x) which suffice for theorems similar 
to Theorem 6.1 may be found in [1] and [6]. 

When A(t) is not constant in (6.1), then Theorem 6.1 need no 
longer necessarily hold. However, if the linear portion x' = A(t) x 
of (6.1) is uniformly stable, then we have the following theorem. 

Theorem 6.2 . If the linear system x' = A(t) x is uniformly stable 
and f(t,x) satisfies 

(6.7i) f(t,x) defined and continuous for t ^ 0, 

I |x ! I <5 and 

(6.7ii) ||f(t,x)|| a(t) ||x||, where a(t) > 0 and 

00 

/ a(t) dt < 

0 

then the solution z(t) e 0 of x' = A(t) x + f(t,x) is uniformly 
stable on [0,«>) . 

Proof: Let X(t) be the fundamental matrix of the linear system 

x' = A(t) X. Then there exists M > 0 such that | |X(t) X"^(s)|| < M 
for t ^ s ^ 0. As in the proof of Theorem 6.1, it can be shown that 
any solution x(t;0,x^) of x’ = A(t) x + f(t,x) exists and is defined 
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on [0,~). Then assuming this, the solution x(t;0,x^) can be written 

1 ^1 

x(t) = X(t)X" (s)x(s) + /_ X(t)X“ (T)f{x,x(T)) dx, 

0 

where s ^ 0 is fixed, so that 

I |x(t)M £ M| |x(s)l 1 + / M a(x)| lx(x)| Idx. 

0 



Applying Gronwall's Lemma yields 



I lx(t) I !_< M] |x(s) 1 1 exn / M a(x) dx 

0 



and thus 

00 

I IxCt) II £ M| |x(s) I I exp / M a(x) dx i K] |x(s)| | 

0 



where K is a convenient constant. The above relation is valid for any 
s >. 0 so that given e>0, ||x(s)|| <■!■ for some s ^ 0 implies 
I |x(t)l 1 < e for all t ^ s so that the trivial solution z(t) = 0 of 
x' = A(t) X + f(t,x) is uniformly stable. 

Corollary 1. Under the same hypothesis as Theorem 6.2, if the 
linear system x’ = A(t) x is in addition asymptotically stable, then 
so is the trivial solution z(t) = 0 of x* = A(t) x + f(t,x). 

Proof: The hypotheses imply that z(t) h 0 is a uniformly stable 

solution of x' = A(t) x + f(t,x). If the linear system is in addition 
asymptotically stable, then ||X(t)|l -> 0 as t->« where X(t) is the 
fundamental matrix of the linear system. Therefore, if x(t;0,x^) is any 
solution to x' = A(t) x + f(t,x), then given e>0 there exists T such 
that 1 |X(t)l 1 < ^ fort^T. Then from 



we have 



t 

x(t) = X(t)x„ + f X(t)X 
° 0 



-1 



(s)f(s,x(s))ds, 



l|x{t)l| < llX(t)M llx l| + ) 

° 0 



M a(s) 1 lx(s) 1 1 ds. 
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where | |X(t) X"^ (s) | | <_ M for t ^ s ^0. Then for t ^ T we obtain 

t 

I |xCt) I I 1 e + M f a(s) 1 |xCs) I I ds 
0 

and applying Gronwall's Lemma, 

t 

I |x(t) 1 1 < e exp / M a(s) ds , 

0 

oo 

so that I |x(t) 1 1 ^ e exp / Ma(s)ds<_eK 

0 

where K is a suitable constant. Since e >0 is arbitrary and K is 
independent of c and T, we conclude that ||x(t)||-*-0 as t-^“ and 
that the trivial solution z(t) ^ 0 of x' = A(t) x + f(t,x) is 
uniformly and asymptotically stable. 

From the proof of Theorem 6.1, we note that when the linear 
system x' = A(t) x is uniformly stable and f(t,x) satisfies 
(6.7i,ii), then all solutions of (6.1) remain bounded. However, 
in the event the linear system is restrictively stable, then the 
solutions of (6.1) are actually bounded by those of the linear 
system as the following theorem illustrates. 

Theorem 6.3 . If the system x' = A(t) x is restrictively stable 
and f(t,x) satisfies (6.7i,ii), then every solution x(t;0,x^) 
of x' = A(t) X + f(t,x) satisfies x(t) = 0(||X(t)||), i.e., 

I |x(t)| I 1 K I |X(t)| I I |X|^| I where K > 0 is constant. 

Proof: Let X(t) be the fundamental matrix of the linear system 

x' = A(t) X. Then there exists M > 0 such that ||X(t)|| and 
I |X”^('t)l I ^ ^ t ^ 0. Let x(t;0,x^) be a solution of 
x' = A(t) X + f(t,x). Then from 

t , 

x(t) = X(t)x^ + X(t) X"'(s) f(s,x(s)) ds. 
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we obtain 



x(t) 

xity 




t 

+ M / a(s) 
0 



I |xCs) I I ds 




+ / a(s) 

0 



x(s) 

x(?y 



ds. 



Applying Gronwall's Lemma yields 



x(t) 

KtT 



-''O' 



exn / 
0 



a(s) ds 



< 1 |x II exp / M a(s) ds < K| |x | | 

-‘'o Q - ''o''* 

where K is a convenient constant. Thus M x(t))ll IK llx|l Mx(t) 1 1 
which completes the proof. 

Theorem 6.3 shows that when the linear system is restrictively 
stable, the contribuition of the non-linear term f(t,x) to the solutions 
of (6.1) is small. Golomb [4] has given conditions under which bounded- 
ness and stability may be determined in systems of the form in which 
f(t,x) satisfies 

(6.8i) f(t,x) is continuous in x a.e. for t ^0 

and is integrable in t for fixed x and 

(6.8ii ) I |f (t,x) II 1 a(t) u( | [x | | ) for t ^ 0, where 

a(t) is positive and integrable and u(r) 
defined for r> 0 is a positive non-decreasing 
continuous function which is submultiplica- 
tive, i.e., 
u(r s) ^ u(r) u(s) , 
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without requirinq uniform stability in the linear system, essentially 
qeneralizinq Theorems 6.2 and 6.3. Other evaluations of solutions to 
(6.1) may be found in the same paper. 
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VII. CONCLUSIONS 



We have considered stability in the sense of Lyapunov for 
systems of ordinary differential equations with oarticular emohasis 
on linear systems. The theory for these systems is quite advanced 
today and in Sections III, IV, and V we have considered several of 
the more important results. Though by no means exhaustive, much of 
the material in these sections comes from varied sources and it is 
felt that the presentation there offers a more complete treatment 
than is usually found elsewhere. In addition, a new set of sufficient 
conditions for a stability comparison theorem have been formulated 
which to the best of the author's knowledge have not aoneared else- 
where. Though the treatment in Section VI is somewhat brief, it is 
intended for the most part to complement the rest of the paper by 
illustrating some of the more useful properties exhibited by uniformly 
and restrictively stable systems. The theory of stability for non- 
linear non-autonomous systems is still quite incomplete and is today 
an area of extensive active research. 

Some other aspects of stability theory not considered here are 
the Poincare-Bendixson theory for autonomous systems, the second 
method of Lyapunov, and various analytical approaches such as the 
method of Lindstedt, Krylov and Bogolyubov, and of Van Der Pol. 

As mentioned earlier, the approach to the problem of stability of 
a system of differential equations varies with both the form and 
order of the system. This paper has considered but one of these 
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approaches, one which is particularly useful in the treatment of 
linear systems. For a treatment of some of the other areas in the 
theory of stability of ordinary differential equations, we refer the 
reader to several of the excellent texts on stability, including 
[2], [5] and [7]. 
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